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ABSTRACT 



We consider the general behaviour of cosmologies in Brans-Dicke theory where the dilaton 
is self-interacting via a potential V(<&). We show that the general radiation universe is 
a two-dimensional dynamical system whereas the dust or false vacuum universe is three- 
dimensional. This is in contrast to the non-interacting dilaton which has uniformly a 
two-dimensional phase space. We find the phase spaces in each case and the general 
behaviour of the cosmologies. 
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1. Introduction 

Einstein's theory of General Relativity is extremely successful at describing the dy- 
namics of our solar system, and indeed the observable universe, nonetheless, the realization 
that general relativity probably does not describe gravity accurately at all scales has led 
to various alternatives being explored, most notably the scalar-tensor family of gravity 
theories pioneered by Jordan, Brans and Dicke (JBD) [|IJ, and the gravitational lagrangian 
inspired by low energy string theory, ||. In fact, these two theories are not unrelated, 
since the low energy effective action for bosonic string theory generically takes the form 

S = J d'xV^e-^ (r + 4(V0) 2 - iff^j (1.1) 

where is the dilaton, and H^ v \ is the field strength of the two form B^; comparing this 
with the JBD action 



Sjbd = / d 4 xy/=g 



(Vd>) 
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$R-u± '— + 16tt£ 



(1.2) 



shows that the scalar-tensor sectors are identical if u = — 1. 

The phenomenological importance of such extended gravity theories to modern cos- 
mology was powerfully highlighted by La and Steinhardt ||, who suggested that by using 
JBD theory instead of general relativity, the "graceful exit" problem of old inflationary 
cosmology Q might be ameliorated. It was rapidly realised || that this original scenario 
required a value of the JBD parameter uj which was in conflict with the observational 
limits || . This then led to the development of various alternatives , however, a general 
study of the extended inflationary ethos || indicates that in order to satisfy observational 
constraints rather contrived models are required. Nonetheless, extended inflationary ideas 
still survive in some recent inflationary models (e.g. 0). 

The prominence of string theory as a theory of everything, in particular a quantum 
theory of gravity, means that we must examine its consequences in regimes where it de- 
parts from general relativity, in particular we expect that the early universe might display 
"stringy" qualities; Jl(J gives a selection of articles dealing with precisely this problem. 
Correspondingly, by field redefinitions, one can instead investigate the cosmological impli- 
cations of JBD theories, which several recent authors have done in a variety of ways ||11|1 , 
|EJ, |}I3[] , HUj . Most of the recent studies of the JBD theories have focussed on the quali- 
tative behaviour of cosmologies in pure JBD theory, or JBD plus a cosmological constant^. 



' It should be stressed that the 'cosmological constant' in this case amounts to setting R 
R + A in the JBD action. 
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However, since it is generally believed that the string dilaton must be massive in order 
to anchor it to the current observed Newton's constant, it seems likely that more realistic 
JBD models will include a potential for the dilaton: 

S = S JBD - J d*x^V($) (1.3) 

The specific form of this potential is unknown, but it seems reasonable that at high tem- 
peratures it might look something like Vo$ 2n . In this paper therefore, we are interested 
in the qualitative behaviour of JBD cosmologies with just such a form for the potential, 
focusing in particular on n = 1: V(Q) = Vo$ 2 . (The cosmological constant models in- 
vestigated in |13|] have V(Q) = A$.) Clearly this is not a particularly realistic potential 



for the current time, since it has a minimum at $ = 0, which would correspond to an 
infinitely strong gravitational coupling, (although as we will show V(<&) is not necessarily 
to be viewed as a canonical scalar field potential) however, the model potential ought to 
provide some insight on the cosmological effects of a self-interacting dilaton. 



Our starting point will be the recent papers of Kolitch and others|IJ], who showed 
that the field equations for cosmological models could be reduced to a two-dimensional 
dynamical system for any reasonable perfect fluid matter source in pure JBD theory, 
and also for spatially flat models with a non-zero cosmological constant. In the presence 
of the scalar potential, the cosmological dynamical system turns out to be generically 
three-dimensional, apart from radiation dominated universes in which the system once 
more reduces to two-dimensions. The layout of the paper is as follows. In the next 
section we analyse the cosmological equations of motion following from ( |1.3| ) and derive 
the corresponding minimal dynamical system for a spherically symmetric perfect fluid 
cosmology. In the third section we consider the example of empty universes as a means 
of testing the method and exploring the effect of the potential on the purely gravitational 
theory, we also consider radiation universes, since these too turn out to have a particularly 
simple form qualitatively similar to the vacuum spacetimes. We also consider conformally 
transforming to the so-called Einstein frame to motivate some of these features. In section 
four we present an analysis of dust and false vacuum models. Finally, we sum up our 
results and conclude in section five. We use a mostly plus signature for the metric. 



2. The cosmological equations. 

In this section we present a brief derivation of the dynamical system representing a 
spherically symmetric cosmology in JBD theory with the scalar potential. Varying the 
action ( |1.3| ) with respect to the metric and scalar field yields the following equations of 
motion: 



UJ 
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Contracting ( pTl| ) and substituting in (|2.2| ) gives 



dV($) 



8tt m 1 
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(2.1) 
(2.2) 



(2.3) 



as a reduced equation of motion for $. 

Following the cosmological principle, we will assume a standard Friedman-Robertson- 
Walker (FRW) form for the metric: 



ds' 



-dt 2 + a 2 (t) 



dr A 



1 — kr 2 



+ r 2 (d6> 2 + sin 2 6d(j) 2 ) 



(2.4) 



Substituting into (|2.1|) and (|2.3|) gives the following equations 
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for the universe. Note that we have assumed the stress energy takes the form of a perfect 
fluid, T ao =diag(p, p, p, p) in an orthonormal frame. This satisfies the conservation equation 



a 



p = _3_(p + p ) 



a 



(2.6) 



Following jnj we transform to conformal time 

dt 



n 



a(t) 



(2.7) 



and we take an equation of state 

p=( 7 -l)p (2.8) 

for the fluid. Note that this implies 

p oc a" 37 (2.9) 



from (12.61 ) 

"i n (T 

drj 



Denoting -4- by a prime, and defining 



X 



2u + 3 1 1/2 _ ^ 

12 ¥ ~ ¥ 



a' 

y = 7 + 2* 



(2.10) 



we see that for the potential V($) = Vo$ 2 , the latter two equations of (|2.5|) become: 

X' = -2XF+ (l-3 7 /4) ^ (2.11) 

y 2 + fc = X 2 + ^ + M^! ( 2 .i 2) 
3$ 6 y 1 

we may then, | 14fl , differentiate this constraint equation for y, using the X-equation and 



the equation of motion for p to finally obtain 



X' = -2XY + (1 " ^ 7/4) (y 2 + fc - x 2 - I) 

A 6 

y = (i _ 3 7 /2)(y 2 + k - x 2 ) - 2X 2 + ^ ( 2 - 13 ) 
z' = 2zy 

where Z = Vo$a 2 , as the three-dimensional dynamical system representing the general 
perfect fluid cosmo logical models for a JBD theory with the scalar potential Vo$ 2 . (The 
method can be generalised for V(<&) = Vo$ n in which case Z = Vo$ n_1 a 2 , and X' acquires 
additional Z-terms.) 

Note that (|2.13|) represents the most general case scenario. If Vq = 0, the Z-equation 
disappears and we recover the Kolitch and Eardley scenario [14]. Also, if p = (as we will 
consider in the next section) ( |2.12j ) can be used to eliminate Z, again leading to a two- 
dimensional dynamical system, and in fact, we can also reduce ( |2.13|) to a two-dimensional 
system in a radiation universe, since p oc a -4 which allows us to eliminate X using ( |2.12|) . 



Before proceeding to an analysis of ( |2.13 ) for the various cosmologies, we conclude 
this section with a few general remarks. First, note that in order for the Ricci term in the 
action (1.2) to have the correct sign, we require $ > 0, i.e. Z > 0. From ( |2.13| ) we note 
that since Z'\z=o = 0, the dynamical system trajectories will not cross Z = 0. We will 
also take u > —3/2, as is conventional, which implies that A > 0; for the case of string 
theory, A = ^K=. For a physical cosmology, we will require p > 0, i.e., 

Y 2 + k-X 2 ~->0 (2.14) 
6 

from the constraint equation ( |2.12|) . Using the dynamical system equations, ( 2.13 ), verifies 
that [Y 2 + k — X 2 — ■§■] =0 and hence trajectories do not cross p = 0. In other 

words, a positive energy trajectory will remain a positive energy trajectory. 

Finally, it will be of interest to note whether the cosmological models corresponding 
to the dynamical system trajectories are expanding or contracting. From ( |2.10| ) we note 
that for expansion Y — X/2A > 0, therefore we will be interested in where this surface 
is relative to the physical regions and the possible trajectories. It will mostly be the case 
that trajectories will lie on one side or the other of this surface, and hence will represent 
cosmologies that are either eternally expanding or contracting, however, particularly for 
small A, it may well be the case that some trajectories cross this line, in which case they 
will correspond to cosmological solutions which 'bounce' in the Brans-Dicke frame, that is, 
they start off contracting, reach a minimum size and re-expand (or vice- versa). However, 
it should be pointed out that universes which bounce in the JBD frame, may well not be 
true bounce universes when viewed from the 'Einstein' frame, a conformally related metric 
in which the gravitational part of the action appears in Einstein form. For a universe to 
bounce in the Einstein frame requires Y to change sign, something which is only possible 
in k = 1 cosmologies. 



3. Vacuum and radiation cosmologies. 

We consider vacuum cosmologies first, since these should display the main features of 
the effect of the potential, however, we also include radiation universes since these too are 
effectively two-dimensional and are exactly soluble. We first derive the two-dimensional 
dynamical system governing the radiation/vacuum cosmologies, determining the form of 
their solution, and plotting the curves on the two-dimensional phase plane. Using also 
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the original X, Y, Z-variables, we identify the various qualitative behaviours possible, re- 
interpreting key solutions in terms of the original cosmo logical parameters a(t). 
Setting 7 = 4/3 in (|2.13|) gives 



X' = -2XY 



Y> = -(Y 2 + k + X 2 ) + 
Z' = 2ZY 



(3.1) 



The key feature that enables radiation (and hence also vacuum) universes to be simplified 
is that for 7 = 4/3, ( |2.9| ) implies p oc a -4 , which means that the constraint equation ( |2.12j ) 
reads 

Y 2 + k 



SirpoVo 2 Z 
~^Z~ +X + 6 



This in turn allows us to decouple X from (|3.1|) leading to 



Y' = -2{Y 2 + k) + 



87tp V Z 



3Z 



+ 



In order to reduce this to a standard form, we write 



W = YZ (=Z'/2) 



and obtain 



Z' = 2W 

, Z 2 B 

W = 2kZ + — 

2 2 

where we have written B = 167rpoVo/3 for convenience. 
It is not difficult to see that 
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(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



satisfies ( p.5|) , where C is a non-negative constant from (|3.2|) . Since W = Z'/2, (|3.6| ) is in 
fact an elliptic equation, and hence Z(rj) can be written in closed form: 

Z3 — Z\ 



Zx + 



6 '/' 



Z2 — Z-1 
Z3 — Z1 



(3.7) 



where the Zi are the roots of the cubic W 2 {Z). However, since this leads to rather involved 
general expressions for $ and a in terms of integrals of elliptic functions, it proves to be 
more illuminating to proceed with the qualitative picture. 
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First we identify the critical points of the (Z, VF)-system, i.e. points at which W = 
Z' = which represent an equilibrium solution of this system of equations. From ( |3.5| ), 
(|3.6|) , these are 

P+ = (2 + V4^B, 0) ; with C = C (B), k = 1 

(3.8) 

P 0> _ = (0, 0) ; B = C = 0, k = 0, -1 

We will return to the classification of these critical points later. 

Next, note that if W 2 (Z) is always positive then W is always defined. This is clearly 
true for all B and C when k = 0,-1, and for B > 3 if k = 1. If, however, W 2 (Z) < at any 
point, then W cannot be defined and the corresponding trajectory will consist of disjoint 
segments. This means that there are three main forms that the (Z, W) phase plane can 
take, which are illustrated in figure 1, according to whether the cubic W 2 (Z) has three, 
two or one real root(s) for the C = trajectory. (Figure Id is the special case B = 
of k = 1.) Although these figures are shown for k = 1, the diagrams are qualitatively 
the same for k = 0,-1 with these cases generically taking the qualitative form of figure 
lc, although for B — (vacuum) the figure includes a critical point at the origin. Not 
surprisingly the three different pictures translate into three (slightly) different qualitative 
behaviours for the cosmologies. 

Before turning to this however, we first investigate a little further the parameter ranges 
for each option. Note that W 2 (Z) always has one non-positive root, Z\ < 0. Whether or 
not it has additional real roots depends on whether the discriminant 

D = B 3 -3B 2 k 2 + 18kBC + 9C 2 -48k 3 C (3.9) 

is negative or zero. 

• If D < for C = 0, then D < for some range of C < Co, say, where Co is the critical 
value of C for which D(B, C ) = 0: 

Co = |fe 3 - kB + \(Ak 2 -B) 3/2 (3.10) 

Therefore, for C < Co, W 2 (Z) has two additional positive roots for some range of C and 
we have diagram la. This will be the case for k = 1, B < 3,C < Cq. A graph of Cq(B) is 
shown in figure 2. 
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Fig. 1: The (Z,W) phase plane for the radiation cosmologies. Critical points are 
indicated by a disc, the disallowed C < and Z < regions are shaded out. The 
plots are shown for k = 1 cosmologies, but the general form of k — 0, — 1 cosmologies 
is qualitatively the same as figure lc. 
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Figure 2: k = 1, D < 




Fig. 2: A plot of the parameter region in (B, C) space in which the discriminant, 
D, is negative. When D can be negative the (Z, W) phase plane takes the form of 
figure la and we have bounce cosmologies. This can only happen for k = 1. 



If D = for C = then W 2 (Z) has one (repeated) root for C = 0. If k = 1, B = then 

9 



D is negative on (0, Co) and we have diagram Id, otherwise all other trajectories have no 
roots and are continuous. Examining fl3.9|) shows that B = 3k 2 , for D = C = 0, and the 
repeated root is therefore 

Z + = 3k ; Z Q = (3.11) 

according to whether B = 3k 2 or B = 0. If k = 1, B = 3, Z + = 3, and we have diagram 
lb. If k = 0, —1, we can only have B = 0, and we get diagram lc (with the critical point 
at the origin as already mentioned). 

• If D > for C = then W 2 (Z) has no zeros and we have diagram lc. This will be the 
case for all positive B for k = 0, — 1, and for B > 3 if k = 1. 

In terms of interpreting cosmological solutions it is easier to examine the original X, Y 
variables. Using ( |3.2| ), (|3.4| ), and ( |3.6| ) it is easy to see that 

X 2 - ° 

v2 J „ B C < 312 > 
6 + 2Z + Z? 

which can be used to produce parametric plots of the (X, Y) plane. Although these are 
broadly similar for radiation and vacuum spacetimes, the vacuum plot is a true two- 
dimensional dynamical system whereas once the energy density of radiation is nonzero, 
the plot becomes a projection and contains apparent crossings of trajectories. 

In translating from the (Z, W) plane to the (X, Y) plane there will be, broadly speak- 
ing, three qualitatively different pictures. 

• Figure la contains trajectories along which D < 0, and hence Y has a zero. If we order 
the roots of W 2 (Z) as Z\ < Z^ < Z3, then these trajectories have two branches: the 
Z G [0, Z2] branch and Z > Z3 branch. For Z > Z3 the trajectory starts at (0, —00) with 



Z infinite, and crosses the axis at X = ±v CA/Z3 where Z is at a minimum, approaching 
(0, 00), Z — > 00. Since X always has the same sign, $ is either monotonically increasing or 
decreasing along such trajectories. This means that it must be a (77) that is causing Z to be 
infinite, reducing to a minimum, and increasing again, i.e. this is a bounce cosmology both 
in the Einstein and the JBD frame. For the other branch, Z starts at zero, is maximized at 
Z2, returning to zero again. The corresponding trajectories start at (±00, 00) and tend to 
(±00, —00). In this case, X once again has a fixed sign, however, it is no longer small, hence 
$ either increases from to 00 or vice versa. These solutions will also be bounce universes 
in either frame. For the critical D = trajectory there are again two branches, Z > Zq 
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and Z < Zq, where Zq = Z2 = Z^ is the repeated root. Each of these branches corresponds 
to four trajectories, each of which interpolates between the images of the critical point P+ 
in the (X, Y) plane : ( ± ^^, 0), and (0, ±00) for Z > Z , or (±00, ±00) for Z < Z . For 
D > there are no zeros of Y and trajectories interpolate between {(±00, ±00), Z = 0}, 
and {(0, ±00), Z = 00}. 

• Figure lb has one D = trajectory which terminates at (0,0); all others are of the D > 
form. 

• Figure lc has D > for all trajectories hence no trajectories cross or terminate on the 
X-axis and the phase plane consists of two disjoint parts with a 'barrier' across the plane 
separating positive and negative Y. 

Such are the broad features of the phase planes, however, since vacuum cosmologies 
are a special case, we analyse them first before proceeding to the general case. 

3.1. The vacuum cosmologies. 

Although we have solutions for X and Y, we first show how they form a two- 
dimensional dynamical system. For p = 0, ( |2.12| ) can be used to eliminate Z from ( |2.13| ) 
leaving 

X' = -2XY 

Y = Y + k — 3X 2 

independent of the value of A. This dynamical system has four possible critical points: 



5 0ll = (0,±V=fe) <- P_,P ; 



(3.14) 



S 2 ,3=(±\/g,0) <- P+,P 



depending on k. These points are classified in table 1, and the phase planes are shown in 
figure 3. 



TABLE 1: The critical points of a vacuum cosmology. 




k= -1 


k = 


k = +l 


£0,1 


(0,±1), SADDLE 


(0,0), DEGENERATE 




£2,3 




(0,0), DEGENERATE 


a/5,o) ' SADDLE 
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Figure 3: The vacuum cosmologies 
Figure 3a: k = Figure 3b: k = 1 




Fig. 3: The (X, Y) phase plane diagrams for the vacuum cosmologies. The shaded 
areas represent the disallowed regions where Z < 0. Critical points are indicated 
with a dot. Figure 3d shows a closeup of the critical points in figure 3b, the grey 
line represents a' = for A = 0.75. 
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In order to interpret the phase diagrams, we must translate from the X, Y variables 
to a cosmological solution $(£),a(£). We begin with the critical points. 

(i) Critical point solutions. 

• k = 0: Critical point (0,0): For X = Y = 0, a is constant and we observe from ( |2.12| ) 
that $ = 0. Hence this "cosmology" is simply Minkowski spacetime. 

• k = 1: Critical points (^|?0): Here the constraint equation ( |2.12|) gives Z = 4, a 
constant. Solving for a then gives 



a(t) = 

4lf ( 3 - 15) 

m = w 

This now has a non-trivial JBD scalar field. The upper (lower) branch corresponds to the 
+(— ) critical point, has t G (— oo,0) ( t G (0, oo)) and is collapsing to (expanding from) a 
singularity. 



• k = — 1: Critical points (0, ±1): Here (|2.12| ) gives $ = again, however, a' /a = ±1 is 
non-zero. Integration gives a oc e ±r] hence a = \t\. The +(— ) root corresponds to t > (<)0 
and corresponds to an expanding (contracting) universe, however it should be pointed out 
that this is a coordinate transform of Minkowski spacetime, and is known as the Milne 
universe. 

In addition to the critical point solutions, there are also some simple solutions to 

(EH). 

(ii) Constant scalar exact solutions. 

An examination of ( |3.13| ) readily reveals that the l"-axis (or segments thereof) repre- 
sents a solution. Inputting X = gives 

{tan r\ k = 1 
f k = (3.16) 

— coth?] k = — 1 

• k = 0: The y-axis consists of two solution segments, Y > and Y < 0, for which 

a(t) = e^ ; $ = ^ (3.17) 

• k = 1: The whole of the y-axis is a trajectory, which corresponds to the cosmological 
solution 

6H 2 

a(t) = H' 1 cosh Ht ; $ = —- (3.18) 

Vq 
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• k = — 1: There are again two solution segments, since (|2.12 ) requires Y 2 > 1. The 
cosmological solutions are given by 

ctt2 

a{t) = H' 1 sinh H\t\ ; $ = — - (3.19) 

Vq 

Note that these are the de-Sitter family of solutions. 
(Hi) Attractor solutions. 

Finally, it is apparent from figure 1 that there are three main attractors: 
A = {X = 0, Y -> oo} 
A± = {-Y oc ±X -> oo} 

We have already discussed the attractor Ao which corresponds to the de-Sitter family 
of universes, therefore it only remains to discuss A±. Clearly, from ( [3.1 2|) Z — > along 
each of these trajectories, and to leading order we find: 

Y oc — , Z oc — n , with rj — > 0~ 
2r) 

a(t)ac\t\V&, *oc|t|^, with ( * "> 0_ for 6 ^ ± 1 > °" ( 3 - 20 ) 

w 11 11 1 1 -> oo for 6A - 1< 0. 

a(t) oc e Ht , $ oc e" 3Ht , with t -> oo for 6A = 1 

Thus the attractor always corresponds to a (power law) contracting universe with a 
diverging JBD scalar. The A- attractor has a vanishing JBD scalar and is contracting 
for A > 1/2, expanding for A < 1/2. This behaviour can be seen clearly from considering 
the a' = line Y = X/2A. For A > 1/2, this line lies in the physically inaccessible p < 
regions for k = 0,-1, only briefly crossing the p > region for k = 1 giving rise to a 
small number of bounce solutions in that case. For A < 1/2, the a' = line lies (almost) 
entirely in the allowed p > region, with a small exception now for k = — 1 where it 
crosses the p < barrier. This means that there are many bounce solutions for A < 1/2. 
In particular, the vacuum attractor A- now becomes a late time expanding solution. 

We are now in a position to state the general behaviour of a vacuum cosmology. Apart 
from a few special k = 1 trajectories which cross the X-axis, the late time behaviour of a 
vacuum universe is determined by which quadrant of the (X, Y) plane it lies in. Trajectories 
with Y > attract to Aq - the de-Sitter universes, and those with Y < attract to A± 
according to whether X = ±|X|. The exceptions to this rule are a small set of k = 1 
cosmologies which correspond to trajectories crossing the X-axis. These correspond to 
D < trajectories and either have Z > Z3 or < Z < Z2. The former trajectories lie 
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very close to the y-axis and appear to be perturbations of de-Sitter universes with a very 
nearly constant dilaton. The latter trajectories are perturbations of the A± attractors 
either having a very small dilaton (A-) or a very large dilaton (A+). We should also point 
out that in addition to these bouncing cosmologies, for k = 1 there are also a handful of 
vacillating and coasting cosmologies as indicated in the closeup figure Id. A vacillating 
cosmology has two zeros of a' and corresponds to an expanding universe which slows down, 
recontracts for a short period, then re-expands (or vice versa). The coasting cosmology is 
tangent to the a' — line, and represents a universe which expands, slows down and halts 
expansion, then re-expands again (and vice versa). 

3.2. The radiation cosmologies. 

The classification of the radiation cosmologies is extremely similar to the vacuum 
cosmologies, as can be seen from the (X, Y) phase plots shown in figure 4, however, there 
are a few subtle differences that should be highlighted. 

As we have already mentioned, the radiation (X, Y) plot is a projection, and therefore 
contains several features that are a consequence of this. Most obvious is the apparent 
crossing of trajectories. Recall that different trajectories correspond to different choices of 
C in (|3.12[ ), therefore, when Z\ = J^rZ?, the value of X at Z\ and Z2 is the same. If in 



addition Z\ = 6B/Z2, then the Y values are the same. Since Z ranges from zero to infinity 
on those trajectories which have a minimum non-zero value of \Y\, this means that any 
one such trajectory will intersect any other such trajectory once. One implication of this is 
that the boundary of the physically allowed region is no longer a trajectory (or projective 
trajectory) as it was in the vacuum case. The physically allowed region is given from ( |3.2|) : 

Y 2 + k-X 2 = - + — (3.21) 
6 2Z K 1 



In vacuum, Z > gave Y 2 + k > X 2 , hence the shaded regions in figure 3. If, however, 

B_ 

3 



B 7^ 0, the right hand side of Q3.21|) now has a minimum value J| for Z = V3B, i.e. 



/ B 

Y 2 + k-X 2 >\— (3.22) 
V 3 

as shown in figure 4. Therefore, unless we are at the critical point Pq for B = 3, Z cannot 
remain at V3B along a trajectory. This means that each trajectory reaches the minimum 
value at most once, and for Z e (0, 00) trajectories, at least once. Thus the curve given 
by equality in ( |3.22|) forms an envelope for the physical trajectories. 
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Figure 4: The radiation cosmologies 



Figure 4a: k = 1, B = 2 Figure 4b: k = 1, B = 3 




Fig. 4: The (X, Y) phase plane diagrams for the radiation cosmologies. The 
shaded areas represent the disallowed regions where Z < 0. These figures are for 
k = 1, the general figure for k = 0, — 1 taking the appearance of figure 4c. Figure 
4d shows the trajectories for different spatial curvatures with B chosen to give the 
same disallowed region; the solid trajectories are those with k = 1, the dashed, 
k = 0, and the dotted, k = — 1. 
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Another straightforward difference is that for k = 0, —1, there can be no critical 
points, since (|373) implies \Y\, Z > strictly, and hence \Z'\ is strictly positive. Thus 



the intersection of the forbidden region with the Y-axis at (0, (y)*) does not represent a 
critical point, merely a stationary point on the trajectory Y(Z) = J| + For k = 1, 
we must have B < 3 for a critical point, in which case the points correspond to P+ as 
defined in (|3.ij) for the (Z, IY)-plane, and occur at X = y/Co/Z + , where Z + is as defined 
in ( |3.11|) . The critical points are classified in table 2. 



TABLE 2: The critical points of a radiation cosmology. 




k = -1 


k = 


k = +1 


P+ 






(^^b'°)' S<3 ' SADDLE 


Po,- 


VACUUM ONLY 



Finally, note that the 'attractors' A± are subtly different for the radiation case. Al- 
though to leading order A± and the solutions are the same, to sub-leading order, 



in 



i*i + 



B 



4VC 



(3.23) 



Therefore, qualitatively the overall behaviour of the radiation cosmologies is the same 
as the vacuum cosmologies. Universes starting at Y > generically attract to de-Sitter 
universes at late times, and those with Y < to the A± solutions. For k = 1 and 
sufficiently low radiation density (B < 3) there will be a small number of trajectories 
crossing the axis, and hence universes which are a perturbation of the full de-Sitter/ A± 
universes, however, as we consider initial values \Y\ — > oo, these form a smaller and smaller 
set, tending to measure zero. In addition, there will also be a small number of coasting 
and vacillating cosmologies. 

3.3. The Einstein frame. 

Finally, we would like to remark on the naturalness of the de-Sitter solutions within 
this family by making a conformal transformation to the Einstein frame. To do this, we 
set 

9ab = Q^gab (3-24) 
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so that 




(3.25) 



and the gravitational part of the action appears in Einstein form. 

Under such a transformation, the radiation energy-momentum tensor is conformally 
invariant, and the potential transforms as 



i.e. the potential appears as a cosmological constant in the Einstein frame. Thus, for $ 
constant, we would expect de-Sitter solutions. 

4. Dust and false vacuum cosmologies. 



and false vacuum (7 = 0). As before, we identify the critical points of the dynamical 
system, classifying them, and identify the physical regions of parameter space finding 
sample trajectories by numerical integration. Finally, we re-interpret these solutions as 
cosmologies, in terms of the original cosmological parameters <&(£) and a(t). 

4.1. The dust cosmologies. 

For dust, 7=1, and the dynamical system is 




(3.26) 



In this section we analyse the dynamical system ( |2.13| ) for the cases of dust (7 = 1) 




Y' = -l(Y 2 + k) 



3X 2 



Z 




2 



Z' 



2ZY 
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Figure 5: DUST, k=0, A=1 




TABLE 3: The critical points of a dust cosmology. 




k = -1 


k = 


k = +1 


Qo,i 


(0,±1,0) 
SADDLE 


(0,0,0) 
DEGENERATE 




Q 2 ,3 




AS ABOVE 


SADDLE 


<?4,5 


NOT PHYSICAL 


AS ABOVE 





This apparently has six possible critical points: 
Qo,i = (0,±V=M) <- S-0,1 

Q 2 ,3 = (±yf,0,4fc) «-> S 2 , 3 (4.2) 

however, substitution of Q^s into (|2.12| ) indicates that it is a non-physical critical point, 
as it corresponds to p < 0. For /c = 0, these are all coincident at the origin; for k = — 1, 
only Qo,i are relevant, and for k = 1, only Q2,3- Table 3 summarizes the information 
concerning the critical points for the dust filled universe. 

Since the only physical critical points correspond to the vacuum critical points So- 3, 
we can refer to section 3.1 for the solutions to which these correspond, namely, Minkowski 
spacetime or the Milne universe for So,i and ( |3.15| ) for £2,3- 

Clearly, since the dynamical system is now inherently three-dimensional, a two- 
dimensional plot will not convey all the information, however, Since the sign of Z' is 
determined by the sign of Y~, most of the useful information can be gleaned from an (X : Y) 
plot. For the purposes of visualisation, we include fig. 5 which is a three-dimensional 
plot of the A = 1, k = dust system. Note how the trajectories strongly attract to 
Z — > 00, X — > in the Y > region. This corresponds to the de-Sitter attractor Aq of the 
vacuum system. It is not difficult to understand this, since for small X, $ is very nearly 
constant, and a — ► 00 since Z — > 00. In other words, p oc a -3 — > and we rapidly approach 
a vacuum situation. 

For Y < 0, Z — > rapidly along trajectories, and therefore we do not necessarily 
asymptote a vacuum situation. Indeed, since Z ~ 0, we expect, and observe, that the 
cosmologies will take the form of the Kolitch and Eardley dust cosmologies, whose two- 
dimensional dynamical system is obtained from ours by setting Z = 0. 
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Figure 6: The dust cosmologies - a 2-dimensional projection 



Figure 6a: k = 0, A = 1 Figure 6b: k = 0, A A 2 = 1/12 




Fig. 6: The (X, y)-projections of the dust cosmologies. The Ao attractor is evident 
in all of these. Figure 6a corresponds to the projection of figure 5. The grey line 
corresponds to a' = in each case. 



Figure 6 shows a selection of two-dimensional projections of the phase space onto 

21 



(X, Y) plane. Note that A+ is weakly attractive where the divergence of $ mitigates the 
divergence of p in ( |2.12| ) to allow the cosmology to more closely mimic the vacuum scenario. 

To summarize, dust cosmologies which have Y > generally attract to de-Sitter cos- 
mologies. If Y < then A+ appears to weakly attract local trajectories, but other wise 
there is no apparent simple classification of the late time behaviour as in the radiation- 
vacuum cases. For k = 1, as with the radiation/vacuum cosmologies, there exist trajecto- 
ries crossing the X-axis. These correspond to bouncing cosmologies and take one of two 
forms. If Z > Z m i n then the trajectory stays close to the K-axis and is a perturbation of 
the de-Sitter universe. If < Z < Z max then the trajectory start out from (— oo, oo, 0) and 
either asymptotes A+ or A-, the bulk of such trajectories taking the former course. 

4-2. The false vacuum cosmologies. 

For false vacuum, 7 = 0, and the dynamical system is 

X' = -2XY + \(Y 2 + k - X 2 - —) 
A 6 

Y' = Y 2 + k - 3X 2 ( 4 - 3 ) 
Z' = 2ZY 

This has possible critical points 

#o,i = (0,±^=M) <-> 5 ,i 

^ 2l 3 = (±y|,0,4fc) <-> S 2 , 3 (4 .4) 

These are classified in table 4. 

As before, we can refer to the vacuum critical point solutions for Rqs, what is less 
clear is that R^^ are also essentially the same cosmological solutions. To see this note 
that Z = requires either $ = or a = 0, ( |2.12j ) then implies A = or oo respectively 
- a somewhat artificial limit. Setting $ = in order to get a meaningful cosmology then 
gives a oc \t\. This is, in each case, a coordinate transform of Minkowski spacetime and is 
equivalent to the Po critical point of the (Z, VF)-plane. 

As with the dust cosmologies, the false vacuum dynamical system is also three- 
dimensional, and as before, we plot the two-dimensional projection onto the (A, y)-plane 
in figure 7. The false vacuum case however, seems to be the opposite of dust, since for 
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TABLE 4: The critical points of a false vacuum cosmology. 




k = -1 


k = 


k = +1 


Ro,i 


(0,±,0) 

SADDLE 


(0,0,0) 
DEGENERATE 




-^2,3 




AS ABOVE 


V o 

SADDLE 


-^4,5 


i^I^^' 1 ' ); ^<Vv / 3 
A G (0,1/2), FOCUS 
A G (1/2,1/VS), SADDLE 


AS ABOVE 


SADDLE 



V > there is no obvious attractor, while for Y < the attractors show up clearly. 
We have checked this behaviour by compactifying to the Poincare sphere and analysing 
the critical points at infinity which confirms that there is no attractor for Y > 0. Again, 
this is not difficult to understand, since the only difference between the false vacuum dy- 
namical system ([4 . 3| ) and the vacuum one (|3.13| ) is the presence of the extra term in the X' 
equation, which is proportional to Aa 2 /<E>. But (In ^-)' oc (Y — X/A) hence if Y < X/A, 
as it is for both A± for A > 1 or for A+ for all A, then this term is rapidly damped and 
we approach the vacuum situation. For Y > however, this extra term in X' grows and 
we depart increasingly from vacuum at late times. Thus, in distinction to all the other 
cosmologies we have considered, the expanding false vacuum cosmologies do not attract to 
de-Sitter universes at late times. 

5. Conclusions 

In this paper we have presented a dynamical systems analysis of the FRW cosmologies 
in Brans-Dicke theory with a potential for the scalar JBD field. Our analysis was based on 
that of Kolitch and others, |14j], who found that cosmology in pure JBD theory could be 
expressed as a two-dimensional dynamical system. In the presence of a scalar potential, we 
found that unless the matter source was radiation, the dynamical system was inherently 
three-dimensional. In spite of this we were able to see clearly that the asymptotic late 
time behaviour of the expanding cosmologies was generically that of a de-Sitter universe, 
the specifics depending on the topology of the spatial sections of the cosmology, and the 
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Figure 7: The false vacuum cosmologies - a 2-dimensional projection 
Figure 7a: k = 0, A = 1 Figure 7b: k = 0, A A 2 = 1/12 




Figure 7c: k = -1, A = 1 




Figure 7e: k = 1, A = 1 





Figure 7d: k = -1, A A 2 = 1/12 




Figure 7f:k = 1,A A 2 = 1/12 




Fig. 7: The (X, Y)-projections for the false vacuum cosmologies. Note the change 
in the attractors for Y < for the different values of A. The grey line corresponds 
to a' = in each case. 



only exception being that of a false-vacuum source. Since the potential chosen is actually 
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equivalent to an Einstein cosmological constant, these results are in general agreement with 
the statement that general relativity is an attractor for Brans-Dicke theory |12|| , although 
this statement comes with certain caveats, as the false vacuum case illustrates. 

In order to get an idea of when general relativity is a late time attractor, consider the 
intermediate equation for X' , ( |2.11|) , substituting the behaviour of p from (|2.9[) : 

X' = -2XY + (1 - 3 7 /4)^a 2 " 3 ^- 1 (5.1) 

Now, 

(lnfa 2 " 3 "^- 1 ])' oc (2 - 3 7 )F - 3(1 - -y)X/A 

therefore, if 7 > 2/3, [a 2_37 $ _1 ] will decrease along a positive Y trajectory, and X will 
be driven to zero. $ will therefore become constant and we will asymptote a de-Sitter 
universe. For 7 < 2/3 however, we do not expect a late time de-Sitter solution. Since 
p + 3p = (37 — 2)p, 7 < 2/3 corresponds to a matter source which does not satisfy the 
strong energy condition, therefore it appears that general relativity is an attractor for 
matter obeying the SEC. 

Clearly this model is somewhat artificial since it neither anchors $ at a specific value, 
nor behaves as we might expect a scalar potential to behave, however, it does illustrate 
differences that crop up between pure JBD theory and those with a scalar potential. In 
particular the appearance of a strong attractor for expanding universes. However the 
model also displays many similarities with the pure JBD case, for instance, the presence 
of vacillating or coasting cosmologies for k = 1. 

Perhaps the most important conclusion to draw from this work is that for "ordinary 
matter" (i.e. matter obeying the strong energy condition) general relativity is an attractor, 
therefore, independent of the specifics of a dilaton potential, i.e. how and at what energy 
scale the dilaton acquires a mass and gravity becomes Einstein in nature, the universe 
should exit the stringy era of gravity in a reasonably familiar Einstein form. On the other 
hand, at higher energies, when an inflationary phase may have occurred, there is no strong 
attractor (s) - very much in contrast to pure JBD |14|] - hence (as was hinted in ||15|p the 
behaviour of inflationary universes in this case could be strongly affected by the behaviour 
of the dilaton potential. 
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